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Abstract. 

The creation and propagation of sound waves in two-component Bose-Einstein 
condensates (BEC) are investigated and a new method of wave generation in binary 
BEC mixtures is proposed. The method is based on a fast change of the inter- 
species interaction constant and is illustrated for two experimental settings: a drop- like 
condensate immersed into a second large repulsive condensate, and a binary mixture 
of two homogeneous repulsive BEC's. A mathematical model based on the linearized 
coupled Gross-Pitaevskii equations is developed and explicit formulae for the space 
and time dependence of sound waves are provided. Comparison of the analytical and 
numerical results shows excellent agreement, confirming the validity of the proposed 
approach. 
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1. Introduction 

Wave phenomena in Bose-Einstein condensates (BEC) represent a very active field of 
research since the first experimental realization of BEC. Matter wave solitons, shock 
waves and different wave instabilities are among the hot topics of the field presently 
under investigation (for a recent review see pQ). In particular, dispersive shock waves 
in BECs consisting of modulated nonlinear periodic waves which decompose into dark 
soliton trains (e.g., when an initial disturbance is represented by a hole in the density 
distribution) or eventually into small amplitude waves (e.g., when an initial disturbance 
is represented by a hump in the density distribution), have been recently investigated. 
Theory of such waves and soliton trains described by the nonlinear Schrodinger (NLS) 
equation was developed in [21 [31 0] and applied to their formation in a one-component 
BEC in El [7]. Shock waves in discrete NLS models were investigated in [8J, [9] while 
experimental observations of shock waves in BEC were reported in [71 [ID] . Although 
these studies were mainly concentrated on single component BEC, binary BEC mixtures, 
both single species (two distinct spin or hyperfine states of the same atom) [11] and 
double species [12], have been created and the existence of several kinds of nonlinear 
excitations in presence of optical lattices has been demonstrated [T3] . 

In spite of the intrinsic nonlinearity of BECs and the relevance of nonlinear 
excitations for applications, small amplitude wave-packets associated with long- 
wavelength modulations of the condensate (matter sound waves) also appear to be 
important. These waves naturally arise in the linear or quasi-linear regimes, e.g. when 
the propagation of small amplitude disturbances on top of the condensate is considered 
or when the inter-particle interactions (nonlinearity) is "artificially" reduced by means 
of Feshbach resonances. They are generated also in a supersonic flow of BEC past an 
obstacle where they form so called "ship wave" patterns located outside the Mach cone 
[El EEl EEl H7J- Small amplitude wave-packets should not be confused with dispersive 
shocks for which larger disturbances or larger nonlinearity are usually required (this 
is particularly true in the multi-dimensional case). From a mathematical point of 
view they can be characterized as waves which, opposite to solitons, entirely belong 
to the continuum spectrum of the linear eigenvalue problem associated with the NLS 
equation. From a physical point of view, these waves are composed of linear modes which 
have a Bogoliubov-like dispersion relation and which, in a pure quantum description, 
correspond to elementary excitations of the system. The spectrum of such excitations 
characterizes the properties of the two-component BEC superfluid [IB] . 

Propagation of sound waves in single-component harmonically trapped BECs was 
experimentally observed in [lOl IT9] in the form of travelling density perturbations. 
Tightly focused far blue-detuned from atomic transitions laser beam applied to the 
condensate acts as expulsive potential and rapidly pushes atoms from the light field 
maximum. Depending on the laser power and pulse duration, different kinds of waves 
are created in the condensate. Smoothly propagating sound waves are created when the 
laser beam waist is narrow compared to the size of the condensate and the light intensity 
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is small enough. Sound waves in two-component BEC's were experimentally observed in 
[20J. While the advantages of experimental techniques for BECs involving laser beams 
are extremely appealing, there is a crucial drawback when it concerns applications to 
two component condensates (in Ref. [20] the slow light technique was used). The reason 
for this is connected to the fact that different atomic species feel different light-induced 
expulsive potentials (due to optical dipole force) and respond to the laser field in different 
manners. This leads to the appearance of interfering waves which complicate very 
much the process under investigation. From this point of view the search of alternative 
methods to generate matter waves in two component BECs is highly desirable. 

The aim of this paper is twofold. From one side we investigate, both analytically 
and numerically, sound waves in two-component BEC's generated by the time evolution 
of a small initial disturbance. From the other side, we present a new method for the 
generation of different kinds of waves in binary BEC mixtures which is based on the fast 
change of the inter-species interaction constant. In view of the recent progress made in 
the manipulation of the inter-component interactions in binary BEC mixtures [21] by 
means of Feshbach resonances, the proposed method appears to be experimentally quite 
feasible and has the advantage that it relies entirely on the internal physical mechanism 
for the creation of the density disturbance. This method is used to investigate sound 
waves in two experimental settings: i) the case of a drop-like condensate immersed 
into a second large repulsive condensate and ii) the case of a homogeneous mixture of 
two repulsive BEC's both with finite backgrounds. In case i) we first create an initial 
disturbance in the form of a gray or a bright hump (depending on the sign of the 
interaction) on a finite background by adiabatically turning on the interaction between 
the drop-like condensate and the surrounding condensate, and then we rapidly remove 
the inter-species interaction. The energy is then entirely released in the form of sound 
waves in the background condensate. In case ii) the initial disturbance is created by a 
change of the intra- or inter-species interaction constant via optically induced Feshbach 
resonance. To illustrate our approach we numerically simulate the creation and evolution 
of matter sound waves in a two component BEC. A mathematical model based on the 
linearized Gross-Pitaevskii equation is has been developed and explicit formulae for 
the space and time dependence of the sound waves are obtained. Comparison of the 
analytical and numerical results shows a very good agreement, confirming the validity 
of the proposed approach. The obtained results can be of interest for the physics of 
interacting two-component superfluid systems. 

The outline of the paper is as follows. In Sec. II we illustrate the proposed 
approach for the one dimensional setting, by formulating the mathematical model and 
performing numerical simulations. In this case one component is a bright soliton or 
drop-like condensate immersed in the second large repulsive condensate with a finite 
background. Since the first component is on the zero background, after the interspecies 
interaction is turned off, we consider the waves generated in the second (background) 
condensate. Analytical formulae describing the space-time patterns of sound waves in 
coupled BECs are explicitly derived. In Sec. Ill the proposed method is extended to the 
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two dimensional case. In Sec. IV we consider the most interesting and experimentally 
relevant case of sound waves of two component repulsive condensates both with a 
finite background. In this case the inter-species interaction is present during the whole 
evolution of waves in both condensates, exhibiting interference effects in the dynamics. 
In Sec. V the main results are summarized. 



2. One-dimensional condensates 



2.1. Decay of a small disturbance of a stationary state to linear wave packets 

Consider a quasi- ID two-component BEC with attractive interaction in one species 
and repulsive one in the other species. In absence of the interspecies interaction the 
stationary ground state corresponds to a bright soliton in the attractive condensate and 
the Thomas-Fermi (TF) distribution in the repulsive one (which can be considered as 
uniform if the TF radius is much greater than the bright soliton's width). Adiabatic 
turning on the interspecies interaction leads to the formation of a new stationary state 
with a hump or a hole in the repulsive BEC, depending on the sign of the interspecies 
interaction (see Fig. [I]). 
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Figure 1. The stationary state of two coupled ID BECs emerged from the attractive 
(left panel), and repulsive (right panel) inter-species interactions. Zero background 
condensate (solid line) is attractive, and finite background condensate (dashed line) is 
repulsive. 



A fast turning off of the interspecies interaction will then create a disturbed 
stationary state in the repulsive condensate. In this section, we shall study evolution 
of such a perturbation. Since the attractive condensate has a zero background density, 
we can consider only waves created in the repulsive finite background condensate using 
the scalar GP equation. 

In this case dynamics of the repulsive BEC is governed by the GP equation 

iifh + ^xx + - g\^\ 2 )^P = 0, (l) 

where 



A* = 9Po = c 2 



(2) 
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is the chemical potential, p is the density of an undisturbed uniform state, and c denotes 
the sound velocity of long wavelength waves propagation along such a uniform state. 
We imply here that the coupling constant g = gm corresponds to the ID geometry of 
the trap with strong radial confinement (u± ^> uo x ) of the BEC. 
Let the initial state be described by a disturbed wave function 

ipo = y/po + 5ip , \5ipo\ < (3) 

For example, 5ipo can describe the hump or the hole created in the density distribution 
by the interaction with another condensate. Evolution of a small disturbance 5ip is 
governed by a linearized equation (TjQ), 

i8ip t + ^8ip xx - c 2 (5^ + 5^) = 0, (4) 
which should be solved with the initial condition 

H\ t =o = S M X )- ( 5 ) 
It is convenient to separate in this equation the real and imaginary parts by introducing 

5ijj = A + iB, (6) 

so that it transforms to the system 

At + \B XX = 0, . . 

B t - \A XX + 2c 2 A = 0. U 

The variable B can be easily excluded and we arrive at a single linear equation 

Au — c 2 A xx + —A xxxx = 0, (8) 

which can be readily solved by the Fourier method. To this end, we notice that linear 
harmonic waves A oc exp[i(kx — uit)) satisfy, as it should be, to the Bogoliubov dispersion 
law, 



k 2 

co = ±co(k), u){k) = kdc 2 + — . (9) 

Hence the general solution of Eq. ([HD has the form 

Wi ( k y(k x -u, m ™ + / W r 2 ( A; ) e <(te+ w (*)t)^ ) (1Q) 

-oo 27r J-oo Z-K 

where W\^{k) are to be determined from the initial conditions: 
A(x,0) = A (x) = JZlWi(k) + W 2 (k)]e ik *§, 
A t (x,0) = A t0 (x) = -i JZo^( k )[ w i( k ) - W 2 (k)]e ik *§ . 

The inverse Fourier transform gives the equations 

A (k) = W x {k) + W 2 (k), A t0 (k) = -iuj(k)[Wi(k) - W 2 (k)\, 
which readily provide 

Wl {k) = \U Q (k) + ^ ) A tQ (k)Y 

W 2 (k) = \ Uoik) - ^A t0 (k)\ . 



(12) 
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The first equation allows to express A t0 (k) in terms of B (k), A t0 (k) = (k 2 /2)B (k). 
At last, A(k) and B(k) can be related with disturbances of the density and the phase 
in BEC, ip 



Vpo + d~pexp(i8(j)) 
6p 



VPo + Tm + iy/Poty* so that 



.4 



B 



Hence, Eqs. (TT2T) take the form 



4^ 
1 



W2(k) = ^[S PQ (k)-i^SMk) 



(13) 



(14) 



Let us consider, for the sake of definiteness, the case of an initial disturbance of 
density only, 5(f) = 0, so that W\(k) = W 2 (k) = 5 po(k) / (4^/p^) and then 

6p(x,t) = 2y/p^A(x,t) = ^JZ D 6 Po( k ) [e^^+e*^^)*)] dk 

= \ J °° o~Po{k) cos(fcx) cos (tk^J~c 



'c 2 + 



k 2 



dk. 



(15) 



The integrals here can be estimated for large t, as usual, by the method of stationary 
phase. For x > 0, only the first exponent gives contribution to the integral, 



where 



Sp(x,t) 



f(k) 



1 f°° 

— / 5p (k)e ltf{k) dk, x>0, 
4vr J_oo 




(16) 



(17) 



The values of k at the points of stationary phase determined by the condition df /dk = 
are equal to (x > ct) 



k = ±k (x,t), ko(x,t) 
and integration yields 



where 



Sp(x,t) 



d 2 f 



5p (k) 



cos 



2nt 



d 2 f 
dk 2 



c 

72 



tf(k 



x y 
~ct) 



4 + 



(-)\ 



ct\ \ctJ 



1/2 



;is) 



(19) 



dk 2 



k (6c 2 



(4c 2 + fc2)3/2- 



(20) 



Here ko is a function of x and t determined by Eq. f|T8|) . In Fig. [2] we compare the exact 
solution ( [TBI) corresponding to a gaussian initial disturbance 



5p (x) 



x 



exp , 5p (k) = exp 



k 2 a 2 



(21) 



with its approximation f|T9l) . and corroborate these two curves with numerical solution 
of the GP equation where we introduce as initial condition the wave profiles shown in 
the left panel of Fig. [TJ As we see, the asymptotic formula (|T9|) is quite accurate even 
for not very large values of t along almost all wave packet where \x\ > ct. 
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Figure 2. Left panel: Time evolution of the density disturbance in the background 
BEC (p(x) = \4>(x)\ 2 ), when the inter-species interaction constant is set to zero 
(512 = 0) at t = 0, according to GP equation with (722 =.9=1- Right panel: 
Comparison of the exact analytic solution (|15[) . asymptotic approximation (|19p. and 
numerical solution of the GP equation for the density disturbance, corresponding to 
t = 4 and c = 2. The initial disturbance has a gaussian form (|2Tj) with a = 0.95. Space 
dependence for analytical solutions is shown for x > at the moment t = 4. Similar 
calculation gives a symmetric distribution for x < 0. 



3. Two-dimensional condensates 



3.1. Decay of small disturbances to wave packets 



This problem can be solved by the same method which was used in ID case with 
replacement d 2 /dx 2 — > A± = d 2 /dx 2 + d 2 /dy 2 . As a result we obtain the general 
solution in the form 

A(r,t) = J Wl (k) e 'C^«*> * + J W 2 (k)e^ r+ ^j^ 2 , (22) 

where Wi^k) are determined again from the initial conditions, which gives 
^( k ) = i^(^o(k)+^0o(k)) 



^(k) = ik(Wk)-i^o(k)). 



(23) 

For cylindrically symmetric initial disturbance of density only (<50o(k) = 0) we get for 
the wave of density Sp(r,t) = 2 v /p^A(r, t) the expressions similar to Eqs. ffl~5l) . 

5p{r,t) = ^ j™5p {k) [ e Kkr cos V -u(k)t) + e i(kr cob V +u(k)t)] kdkdtpj 

N (24) 



J (kr)kdk, 



where 



J (z) = ^ [°° e izcoB ">d<p (25) 

is the Bessel function of the first kind. 

For fcr>l the integral in (|2^|) can be estimated by the method of stationary phase 
with the use of asymptotic expression for the Bessel function, 



J (kr) « W-?-cos(Ar -tt/4). (26) 
w 7rfcr 
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As a result we get 



5p(r,t) 



sin [tf(k )\ 



2-nJrt 



d 2 / 
dk 2 



(27) 



where 



k {r,t) 



/(*) 



c 

72 



L) 2 _ A + L (L 

ctJ ct\ \ct 



--,1/2 



, r > ct, 




(2f 



(29) 



and |d f/dk 2 \ is given by (1201) . 

In order to compare the analytical predictions with numerical solution of the GP 
equation we need to perform a sequence of simulations as in ID case. However, there is 
an essential difference in the 2D case due to the fact that no stable localized solution for 
the drop-like condensate exists for the 2D GP equation without an external potential. 
To this regard we recall that the attractive 2D NLS equation admits a localized solution, 
the so called Townes soliton, which exist only for a critical value of the norm (number 
of atoms) and is unstable against small fluctuations of the critical norm (overcritical 
norms lead to collapse while undercritical ones lead to background decays). To avoid 
this problem, we prepare the drop-like condensate in the ground state of the parabolic 
trap V(x,y) = v (x 2 + y 2 ) with strength vq = 1, using the imaginary-time relaxation 
method [23]. Then we introduce the obtained ground state solution for the drop-like 
condensate tpi(x,y) and constant amplitude solution for the background condensate 
(-02 {x, y) = 2.0) as initial condition in the coupled GP equation with gi 2 = 0, and slowly 
increase the inter-species interaction constant to a finite value g\2 = —0.5 (see Fig. [3] 
left panel). It is remarkable that after the localized coupled state has been created by 
this method, it remains stable even after the parabolic potential is removed, due to 
the stabilizing effect induced by the interspecies coupling (this state exists only due to 
coupling). As a result the initial state shown in the right panel of Fig. [3] is obtained. 
The density waves in the background condensate can be generated by fast turning off 
the interaction constant gi2 — > 0. Decaying density hump induces circular wave in 
the background BEC traveling outward from the center. We have shown in Fig. H] the 
dependence of 5p on r at fixed value of time (t = 4) for the gaussian initial distribution 



<Wr) 



7TGr 



exp 



6po(k) = exp 



k 2 a 2 



4 



(30) 



with a = 1 calculated according to the exact expression ( 1241) . The asymptotic 
approximation (127]) agrees very well with the exact expression ff24|) at x > ct (not 
shown in figure because the corresponding lines overlap). 



Matter sound waves in two- component Bose-Einstein condensates 



9 




Figure 3. Left panel: Formation of a density hump on the background condensate 
ij}2{x,y), shown at y — cross section, when the inter-species interaction constant is 
slowly raised from zero according to g\2 = — 0.5tanh(5i/i e „d), with t enc t = 100. Due 
to inter-species attraction, the drop-like repulsive condensate ipi(x,y) with a norm 
N = 6ir becomes more narrow. Right panel: The final stationary coupled state of the 
two repulsive BEC's with gu = 1, 922 = 1, with attractive inter-species interaction 
constant 912 = —0.5. 
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Figure 4. Left panel: Density plot p(r) = \tp(x,y)\ 2 of the repulsive background 
condensate at t = 4 after turning off the inter-species interaction constant gi2 = 0, 
according to numerical solution of the GP equation. White stripes correspond to higher 
density. Right panel: Comparison between analytical expression (24), for the density 
disturbance in the background condensate at t = 4 with corresponding numerical 
solution of the GP equation (\ip(x,y = 0)| 2 — po), shown for y = cross section. 
Parameters: po = 9, c = 3, a = 1.0. 

4. Two-component condensates 

Let us now consider the relevant case of a homogeneous binary BEC mixture, when 
both components have a finite background. Two-component condensates are described 



Matter sound waves in two-component Bose-Einstein condensates 10 



in our non-dimensional variables by the equations 

i 8 -lt + \Mi + W - (9n \4>i\ 2 + 9i2 \H 2 Mi = o, 

+ I A ^ 2 + [p 2 - (g 12 |^| 2 + g 22 |^ 2 | 2 )]^ 2 = 0, 



(31) 



where 

Aii = gnPio + fi'i2P2o, /J-2 = gnPio + 922P20 (32) 

are chemical potentials of two species with equilibrium uniform densities pio and p 2 o; 
we keep here the general notation for the Laplacian A which will be specified later 
depending on the number of spatial dimensions in the problem under consideration. 

In disturbed condensates we have ipi — \fP10 + <^0i> ^2 = \J~P~2a + <^02, where 
evolution of disturbances is governed by the linearized equations 

i5ip! + |A5Vi - guPio(Sipi + Sip$) - 9l2y/pl0p20(Slp2 + Sip 2 ) = 0, 
iSip 2 + ^A5ip 2 - gi2^piop2o{tyi + Sipl) - g 2 2p2o{5ip2 + #1) = 0. 
As before, we introduce 

5iJ; 1 = Ai + iB 1 , 5ip 2 = A 2 + iB 2 , (34) 

so that after substitution into (13"31 . separation of real and imaginary parts followed by 
exclusion of Bi and B 2 we arrive at the system 

Ai, tt - gupwAAi + \A 2 A 1 - gi2^Pwp2oAA 2 = 0, 
A 2 ,tt - g22p2o^A 2 + \A 2 A 2 - gi2y/Piop2oAA l = 0. 

For plane waves Ax, A2 oc exp[i(kr — ait)] we reproduce the known dispersion laws 



k 2 

u ± = k\lc 2 ± + - } (36) 



where the sound velocities corresponding to infinite wavelengths are equal to 

r / 1 va 

gilPlO + #22P20 ± y (fl'llPlO - g22p20) 2 + ^g\ 2 pWp20 



The ratios A\j A 2 in these two modes of linear waves are given by 

M 2g 12v /p 10 p 20 

A 2 g22P2o - gnPw ± (c+ - C 2 ) ' 
Thus, the general solution of the system (1351) is given by 



+ j W^e^-^^fiy + / Wu(k)e*C^*> , 

A 2 (r, t) = J W^e^-^^y + J W+{k)e 
+ j W 21 (k)e^-- 0, + j W 22 (k)e^^ , 



j(kr+o; + t) 



(37) 



(38) 



(39) 



(40) 



where .D denotes dimension of the space. According to (1381) the Fourier components are 
related with each other by the equations 

Wn 2gi2y / piop2o 



W 21 W 2 % g 2 2p2o - gnPio ± (4 - c 2 ) 



(41) 
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that is only four of them are independent of each other and should be determined from 
the initial conditions. 

To simplify the notation, in what follows we shall consider the case when there exist 
only initial disturbances <5pio(r), 5p 2 o(r) of the densities. Then 

w± = w± = wf, w± = w£ = Wf, (42) 

and these four functions must satisfy the linear system 



A w (k) = 2(W+ + Wr), A 20 (k) = 2(W+ + W 2 ), 



w+ 



2 gi2VPl()P20 



W 1 _ 2g 12y /pu)P20 

922P2()-gilPlO + (c 2 + -ci) ' W 2 ~ 922P20-9llPl0-(c 2 + -C 2 _) ' 



which can be readily solved. Taking into account that 



Aio(k) 



1 



-_S Pl0 (k), A 20 (k) 



2 V / PlO 2y/p~20 

we can find the waves of the densities of the two species. 
In ID case the final formulae read 



<5p2o(k), 



(43) 



(44) 



6pi(x,t) = ± Jjj* ' 5pi{k) coskx 

_|_ ^ _ gilPlO- U22P2U ) 



^ I gllPlQ— P22P20 



Pl °:ff p20 j cos ( tA^c 2 + f 



cos (^tk\J c\ + 



(45) 



cos ( tk\/ci + ^ - cos ( tkJc 2 _ + 4 



k 2 



dk, 



Sp 2 (x,t) = I / °° <5p 2 (A;) cos kx (l - aiiippo ) cos (tky/ ' c\ + f 



+ (l + gllp |Ig 2p20 ) cos UA^c 2 + f 



dA; 



(46) 



r^iWcosb cos ( tkJc\ + f - cos ( tfc Jc 2 + f 



Similar formulae are obtained for two-dimensional cylindrically symmet 
Mr, *) = /o°° 5pi(k)J (kr) (l + cos (tk 

(l - 911P |I| 2P2 ° ) cos (tA^c 2 + 1 



ric waves as: 



(47) 



cos ( tkJc\ + x ) - cos ( t/cv/c 2 + ^ 



kdk ^ 



<W r >*) = ± J °° 5p 2 (k)J (kr 
+ (l + 9iipw~922P2 ^ cog U^^c 2 + f 



gilPl()-322P20 



-^rr^^Jo^r) 



27r(ci-c ; 



) COS + 



(48) 



COS ^tkyjc+ + — COS (^tkyfl 



c 2 +¥ 



It is easy to see that if there is no interaction between two components of BECs, then 
these formulae reduce to the already known results for one-component condensates. 
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Integrals in (I45p - (I48|) can be estimated by the method of stationary phase which 
gives the approximate expressions for disturbances of the densities in the regions x > c + t 
or r > c + t where we can observe the interference pattern of two modes. 

The present setting with homogeneously mixed two repulsive BECs corresponds to 
the most relevant experimental situation. The case of two repulsive ID BECs is akin to 
previously discussed one (see Fig. [1]), if there is a dark soliton in one component. When 
the repulsive inter-species interaction is slowly turned on a density hump emerges in the 
other component, as illustrated in Fig. [51 Linear waves then can be generated by a fast 
turning off the inter-species interaction. 





Figure 5. Formation of a density defect in two coupled ID repulsive BECs with initial 
wave amplitudes ipiix) = 1,i\)%(x) — tanh(x), and gn — 1,<?22 = 1, when the repulsive 
inter-species interaction is slowly turned on according to law gi2 = 0.5 tanh(5t/i en( i), 
with t en d = 100. Dark soliton in one component gives rise to a density hump in the 
other component. 

Another approach for generation of density disturbances in coupled repulsive BECs 
involves local change of the intra- or inter-species interaction constants by a tightly 
focused laser beam via optically induced Feshbach resonance [22] ■ Stationary bright- 
bright, dark-dark and dark-bright coupled localized states on finite backgrounds can 
emerge in mixture BECs when the interaction constants are locally changed. In 
particular, a stationary bright-bright (dark-dark) localized state in all repulsive case 
occurs if the inter-species repulsion is locally reduced (increased). A bright-dark 
localized state emerges when the intra-species repulsion is locally reduced in one of 
the components. 

After the density defect (hump or hole) has been created, a fast change of the inter- 
species interaction coefficient leads to generation of waves in the condensate caused by 
decaying density disturbance. These results are illustrated in Fig. El where comparison 
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between the analytical prediction and numerical simulation of the GPE is also presented. 



0.4 




X X 

Figure 6. Density modulations at t — 8 according to Eqs. (45)-(46) and numerical 
solution of the coupled GPE. Parameters: pio = 1, P20 = 9, a — 1, gn — 1, 522 = 1 5 
.912 = 0.5. 



Since there is no stable localized solution (like dark soliton in ID repulsive case) 
of the GP equation in higher dimensional settings, creation of density disturbance in 
this case also can be done by illuminating the homogeneous mixture BEC by tightly 
focused laser beam designed to change the intra- or inter-species interaction constant 
via optically induced Feshbach resonance. 




-40 -20 



20 40 
X 



40 -20 20 40 
X 



Figure 7. Evolution of waves in a 2D mixture BEC emerged from the decay of a 
Gaussian density disturbance with a = 1 in both components shown at y = cross 
section according to numerical solution of the coupled GP equations. Background 
wave amplitudes of components and nonlinear coefficients are equal to ij)i(x,y) = 1, 
ip2(%,y) = 3, 511 = 1, (J22 = 1, .912 = 0.5. Due to different sound velocities in the 
components radially outward travelling waves quickly separate. 



In Figs. [7H9] the generation, further evolution and pattern of linear waves at 
particular time are depicted for the 2D coupled repulsive BEC's. As can be seen from 
Figs. HI M and [9], the agreement between analytical predictions and numerical solution 
of the GPE is quite good. 



5. Conclusions 



We have investigated sound waves in two-component BECs and proposed a new method 
of wave generation which is based on a fast change of the inter-species interaction 
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Figure 8. Density plots of two components corresponding to the last time section at 
t = 8 of the previous figureJT] In white areas the density is higher. Separation of waves 
in the two components due to different sound velocities is clearly seen. 
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Figure 9. Comparison of density distributions provided by Eqs. (47) and (48) with 
numerical solution of the coupled GP equation for the same parameters as in Fig. [7] 
for t = 8. Relative positions of the wave maximums in the second and first modes 
26/9 ~ 3 well match the ratio of corresponding sound velocities c + /c_ which have the 
same order of magnitude as group velocities of our wave packets built of harmonics 
with wavelengths about healing length. 



constant. To illustrate the idea, we have performed numerical simulations of the creation 
of a density defect (hump or hole) on the background condensate and its evolution into 
matter wavepackets after the fast change of the interspecies interaction constant. This 
method of wave generation was used to investigate both numerically and analytically 
sounds waves arising from the release of a drop-like condensate immersed into a second 
large repulsive condensate and from homogeneous mixture of two repulsive BEC's both 
with finite backgrounds. We demonstrated that propagation of linear waves created by 
this approach is well described by the Gross-Pitaevskii equation. Explicit analytical 
formulae were obtained for the space-time evolution of density waves in coupled BECs, 
starting from the linearized GP equations. The comparison of the analytical results 
with direct numerical simulations of the coupled GPE equations showed an excellent 
agreement. 

Although in this paper we have restricted our considerations to weak density 
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disturbances (e.g. amplitude of the disturbance is small compared to the amplitude of 
the background condensate), the proposed method can be used to investigate also ship 
waves and dispersive shocks in binary BEC mixtures as it will be discussed elsewhere. 
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